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1. Motivations and structure of the paper 

Conformal field theories are a crucial ingredient in both abstract and concrete investigations 
of quantum field theory. They control critical phenomena in condensed matter physics and 
the RG flows of generic quantum field theories [1, 2]. Through AdS/CFT, they even provide 
a framework to study quantum gravity [3]. With some exceptions in two-dimensions, where 
the conformal group is infinite-dimensional, our understanding of generic conformal field 
theories is still quite poor. It is pretty clear that conformal symmetry is a rather restrictive 
constraint on a field theory: the very notion of universality in critical phenomena originates 
from the fact that there are relatively few "simple" conformal field theories which describe 
the infrared behaviour of a large variety of physical systems. The bootstrap program aims 
to use the constraint of conformal symmetry to classify and possibly even solve conformal 
field theories [4, 5]. Recent advances give some hope that the bootstrap strategy could be 
successful even in dimension higher than two, especially when combined with extra input 
from other numerical methods [6, 7]. 

Given the importance of conformal symmetry, it is interesting to consider probes or 
modifications of a theory which preserve a large subgroup of the conformal group. A basic 
example would be a conformal boundary condition, i.e., a boundary condition which is left 
invariant by all conformal transformations which fix the position of the boundary, which 
form an SO{D, 1) subgroup of the S0{D + 1, 1) conformal group of the bulk D-dimensional 
conformal field theory [8, 9]. More generally, we can consider the notion of a conformal 
defect: a d-dimensional defect in a D-dimensional conformal field theory which wraps a 
d-dimensional hyperplane (or a sphere) and is invariant under the SO{d+ 1, 1) x SO{D — d) 
subgroup of the SO{D + 1, 1) conformal group which preserves the hyperplane. Conformal 
defects should play an important role in studying the universal low-energy behaviour of any 



configuration where a quantum field theory is modified or excited in the neighbourhood of 
a large d-dimensional submanifold. It should also be possible to use conformal defects as 
theoretical tools to probe or constrain the properties of conformal field theories. This is 
definitely the case in two-dimensional models, or higher-dimensional superconformal field 
theory, and it may be the case in the context of the bootstrap program as well. See [10] 
for a recent attempt in that direction. 

The purpose of this paper is to study numerically the properties of the twist line defect 
in the critical 3d Ising model, which coincides with the Wilson line defect in the dual Z2 
lattice gauge theory formulation of the model [11] . We will test the hypothesis that the 
twist line defect fiows to a conformal line defect in the continuum limit and determine the 
low-lying spectrum of anomalous dimensions of the operators which live on the defect. Our 
choice of theory and defect is dictated by theoretical and practical considerations. The 3d 
Ising model is currently the basic example of a 3d CFT amenable of a bootstrap analysis. 
The existence and properties of the twist line defects are intimately related to the Z2 fiavor 
symmetry of the Ising model. Thus twist line defects are an example of a conformal defect 
whose existence may encode a crucial property of a CFT. On the practical side, Wilson 
line defects in the Z2 lattice gauge theory are well studied numerically in the confining 
phase as the endpoints for a confining string [12]. This allows us to use well established 
numerical technology. 

The structure of the paper is the following. In section 2 we review the properties of 
conformal defects. In section 3 we define the twist line defects in the Ising model, and 
the basic local operators we will use in correlation functions. In section 4 we present our 
results and conclusions. 

2. Conformal defects 

In the study of quantum field theories, one typically encounters a variety of useful local 
probes and modifications of the underlying theory. The most common example are local 
operators, which probe or modify the theory at a point of space-time. Wilson and 't 
Hooft line operators in gauge theories are classical examples of a probe or modification 
which extends along a line. Boundary conditions or domain walls modify the theory along 
a codimension one locus in space-time. General examples of d-dimensional defects in a 
D-dimensional field theory can be engineered, say, by adding to the Lagrangian of a d- 
dimensional field theory terms which depend on the degrees of freedom of the L'-dimensional 
field theory restricted to the defect . 

If we consider a D-dimensional field theory invariant under the ISO{D) Poincare 
group and fiow to the far infrared, we typically expect to end up with a conformal field 
theory, possibly trivial, topological or free, i.e. a theory which is invariant under the 
SO{D + 1, 1) group of conformal transformations. In a similar fashion, we can consider 
a modification or a massive excitation of the field theory localized near a d-dimensional 
hyperplane and preserving the ISO{d) x SO{D — d) subgroup of the Poincare group which 
fixes the hyperplane. As we flow to the infrared, we do not expect the modification to 
affect the critical properties of the theory in the bulk. We can thus hope that the localized 



modification will either disappear or flow to a "conformal defect", i.e. a defect of the 
confornial field theory, preserving the SO{d+ 1, 1) x SO{D — d) subgroup of the conformal 
group which fixes the hyperplane. 

A priori, the result of the RG fiow is expected to have scale invariance, but full confor- 
mal invariance is a stronger constraint. For unitary quantum field theories, scale invariance 
is expected on general grounds to imply conformal invariance [13]. It is not known if a sim- 
ilar result holds for scale invariant defects in a CFT as well. See [14] for some recent work 
on the subject. Indeed, one basic motivation for the present work was to check conformal 
invariance for a simple example of defect in a non-supersymmetric 3d CFT. 

Although conformal defects are somewhat analogous to local operators, there are some 
important differences. The most obvious difference is that while the set of local opera- 
tors of a field theory is naturally part of the definition of what the theory is, the set of 
higher-dimensional defects which can be inserted in a given conformal field theory can 
be enormous, morally as large as the set of d-dimensional conformal field theories. For 
example, a generic superconformal boundary condition in A^ = 4 4d SYM theory can be 
engineered at weak coupling by gauging a flavor symmetry of a generic three-dimensional 
A^ = 4 SCFT. One may think that the possibility of considering such a variety of confor- 
mal defects is somewhat artificial, and that simple modifications of the theory in the UV 
will lead to a small class of simple defects in the IR. This expectation is incorrect: a very 
simple defect in the UV theory may acquire a very intricate IR dynamics, due to "edge 
excitations" of the bulk theory. Very little is known about possible constraints on how RG 
fiow in the bulk may affect the degrees of freedom living at a defect. 

If the bulk theory is a strongly-coupled CFT in the IR, there is really no well-defined 
separation between degrees of freedom at the defect, and bulk degrees of freedom. The 
closest analogue to studying a "defect conformal field theory" is to look at the set of 
local operators which live at the defect. These local operators have many properties in 
common with operators in a d-dimensional CFT with SO{D — d) fiavor symmetry, with 
one important exception: such a d- dimensional CFT would have a protected stress-tensor 
operator of conformal dimension d. In general, there is no such "defect stress tensor" 
available as a local operator at the defect. 

On the other hand, every conformal defect should support a "displacement operator" , 
which we will denote as D^, i = 1,- ■ ■ ,D — d, which has dimension d + 1 and transforms 
as a vector under the SO{D — d) group of rotations around the defect. Intuitively, the 
displacement operator is something which can be added to the Lagrangian of the theory in 
order to displace the defect in the normal direction, much as the stress tensor is something 
which can be added to the Lagrangian in order to deform the metric of space-time. As 
we can deform the shape of a defect by a local diffeomorphism, there should be a relation 
between the displacement operator and the bulk stress-tensor. 

This relation can be made precise: the displacement operator controls the breaking 
of translation symmetry normal to the defect, and thus enters the stress-tensor Ward 
identities: 

^^T^'i = D'll6{x^) , (2.1) 



where we denote the D-dimensional indices with Greek letters such as ji and the D — d 
transverse indices with latin letters such as i, j, k. This Ward identity makes the protected 
quantum numbers of -D* manifest. Notice that the Ward identity fixes the normalization 
of Z)*, and thus the numerical coefficient in the two-point function 

{D\x)D^m = ^ (2.2) 

is an intrinsic property of the defect. Intuitively, "simple" defects will have a small two- 
point function coefficient. For example, a trivial or topological defect has Cd = 0. 

From the perspective of conformal bootstrap, the correlation functions of the bulk 
CFT can be computed from the knowledge of the spectrum of bulk local operators and of 
the coefficients of three-point functions, whose functional form is determined by conformal 
symmetry. In a similar fashion, correlation functions of defect local operators can be com- 
puted from the knowledge of the spectrum of defect local operators and the coefficients 
of three-point functions of defect local operators. On the other hand, mixed correlation 
functions of local operators in the bulk and on the defect require one extra piece of infor- 
mation: the bulk-to-defect pairing, i.e. the coefficient of two-point functions involving one 
bulk operator and one boundary operator. The functional form of such two-point functions 
is also fixed by conformal invariance: we can use a conformal transformation to send the 
defect operator at infinity, and then use scale transformations, translations and rotations 
to move the bulk operator at whatever location in space-time we want to use as a reference 
point. For example, the correlation function of a scalar bulk local operator O and a scalar 
defect local operator o takes the form 

(O(a;)o(0)) = Cf |x^r2A„|^i|A„-Ao (2.3) 

where |x*| is the distance from the defect, and |x^| the distance from the origin. Indeed, 
when the defect local operator is sent to infinity, the correlation function depends only on 
the transverse distance, as |x*| °~ '^, and then an inversion centered on the origin gives 
the general correlator. 

An alternative, useful point of view is to consider the possible OPE expansions available 
in the system: two bulk operators close to each other can be expanded as a sum of bulk 
local operators sitting at an intermediate location, two defect local operators close to each 
other can be expanded as a sum of defect local operators sitting at an intermediate location 
on the defect, and a single bulk local operator near the defect can be expanded into a sum 
of defect local operators. The latter bulk-to-defect OPE is a good way to use knowledge 
about the bulk local operators to learn about the possible defect local operators. The OPE 
expansion of a bulk operator cannot be empty, otherwise the correlation functions involving 
the defect and that operators would be all zero. For example, if the defect preserves some 
flavor symmetry, there must be a defect local operator for each representation of the flavor 
symmetry for which a bulk local operator exists. 

In this paper we study a conformal defect which belongs to a special class of monodromy 
defects. Monodromy defects can be defined in a CFT which is equipped with a flavor 
symmetry group G. We will focus on the case of discrete G, but most of our considerations 



apply to a continuous flavor group as well. If a conformal field theory has a flavor symmetry, 
we can define a trivial class of topological domain walls T>g associated to elements of the 
flavor group g as follows: a correlation function in the presence of the domain wall is equal 
to the same correlation function without the domain wall, but with all local operators on 
one side of the wall transformed according to g. 

A monodromy defect is defined as any codimension 2 conformal defect on which a T>g 
domain wall can end. Local operators which transform non-trivially under g will be multi- 
valued around the monodromy defect. In particular, this means that the OPE of a bulk 
operator charged under g will contain defect local operators with fractional spin under the 
6*0(2) group of transverse rotations. For example, if (f) is the angular coordinate around 
the defect, r the radial coordinate in the plane perpendicular to the defect, and G = Z2, a 
Z2 odd operator O of conformal dimension A will have OPE 

O(r,0)~j;e^("+i)^^-^o„^i^, (2.4) 

n,a 

involving defect local operators On,a of conformal dimension A^ and half-integral SO{2) 
spin s = n+ 1/2. 

3. A monodromy defect in the Ising model 

Consider the Ising model on a cubic lattice. The Ising model has a Z2 flavor symmetry 
which flips the spin at each site. There is an obvious realization of a Z2 topological 
domain wall in the theory: consider some hypersurface S, and flip the sign of the spin-spin 
interaction for edges which cross S. If S is closed, or extends to infinity, we can simply flip 
all the spins on one side of S, and recover the standard Hamiltonian for the Ising model. 
Similarly, we can deform 5 to a different hypersurface S' , by flipping all the spins in the 
region between S and S' . On the other hand, if S has a boundary, which will generally 
consist of a codimension 2 locus L which does not cross any edges of the lattice, the location 
of L is meaningful, and we obtain a lattice realization of a monodromy defect. Any choice 
of S which is bounded by the same locus L defines the same monodromy defect. 

In the two-dimensional Ising model, the monodromy defect is a local operator, which 
goes in the continuum limit to the disorder operator fi. The OPE between the disorder 
operator ^ and the spin operator a gives a spin 1/2 operator: the free fermion -0 hidden 
in the 2d Ising model. Notice that in the context of the 2d Ising model, the free fermion 
still sits at the end of the topological Z2 domain wall, which in the language of RCFT is 
the topological domain wall labelled by the energy operator e. 

In the three-dimensional Ising model, which is the focus of this paper, the monodromy 
defect is a line operator. The 3d Ising model has a dual description as a Z2 gauge theory, 
in which the monodromy defect is a very fundamental object, i.e. the Wilson loop, and 
has no topological domain wall attached to it. In the Z2 gauge theory description, on the 
other hand, the spin operator is essentially a monopole operator, and sits at the end of a 
topological line defect. The topological line defect acquires a minus sign when crossing the 
Wilson line operator. This implements the anti-periodicity of the spin operator around the 



Wilson loop operator. This is analogue to the behavior of a fundamental Wilson loop and 
a monopole of minimal charge in a 3d gauge theory based on the su{2) Lie algebra. The 
fundamental Wilson loop is allowed if we pick an SU{2) gauge group, the basic monopole is 
allowed if we pick an 5*0(3) gauge group. If we try to include both operators in correlation 
functions, the monopole operator will be anti-periodic around the Wilson loop, and we will 
to place either of the two at the end of some topological defect which keeps track of the 
antiperiodicity. 

After inserting the monodromy line defect in the lattice Ising model, we can tune the 
interaction strength to make the bulk theory critical, and flow to the far infrared. As 
the spin operator is anti-periodic around the defect, the line defect can hardly disappear 
in the IR. It is natural to conjecture that it will flow to a conformal line defect. The 
operators on the monodromy line defect will be labelled by their integral or half-integral 
spin s and conformal dimension A. The displacement operator gives rise to local operators 
D = D^ + iD^ of spin 1, and D = D^ — iD^ of spin (—1), both of dimension 2. 

It is easy to argue that in a theory of a free scalar field, the OPE of the scalar field 
with a monodromy defect would contain spin s = n + 1/2 primary operators of dimension 
|s| + 1/2. Indeed, we can consider a general OPE 

Hx)^Y.fa{z,z)Oa + --- (3.1) 

a 

where z = x^ + ix^ = re"'^ is a complex coordinate in the plane orthogonal to the defect 
(we are taking the defect line along the x^ direction), Oa the primaries on the defect, and 
the ellipsis indicates descendants (derivatives) of the primaries. If we apply the free field 
equation of motion on the OPE and look at the coefficients of primaries, we can ignore 
derivatives along the defect, which give descendants. Thus the OPE coefficients must be 
harmonic functions of ^, z, and the OPE must take the form 

^(X) ~ Y. ^?"+V2. 



^"+^/X+i/2 + c.c +••• (3.2) 



The Ising model is quite close to the theory of a single free scalar field, at least as far 
as conformal dimensions are concerned. Thus we expect the OPE of the spin operator with 
the defect to be dominated by a spin 1/2 operator ip of dimension close to 1, 

a{x) ~ r^i>-^- \e-'^l'^il, + e'-^/^^j + • • • (3.3) 

and that the leading contribution in a spin s = n + 1/2 > sector will be an operator of 
dimension close to n + 1 . 

On the other hand, the OPE of the energy operator e with the defect involves operators 
of even spin, and thus should be dominated by the identity operator, and possibly the 
displacement operator, which we expect to be the operator of lowest dimension in the 
s = 1 sector: 



eix) r^r'^n+r'^-^' 



e~"^D + e"PD 



+ ••• (3.4) 



In a free scalar theory, defect local operators of integral spin could be built as bilinears of 
the scalar field modes, 0n+i/2'^m+i/2^ of dimension equal to |s| + 1. We thus expect the 
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-Monodromy line L 
-Domain wall S 



Figure 1: In our set-up, the domain wall S is the surface across which the links are frustrated. It 
ends on two defect lines. We will mostly consider the proximity of one of such lines, which we take 
to be aligned with the x^ axis. 



Ising model to also include defect local operators of integral spin s and conformal dimension 
close to |s| + 1. It may be possible to understand these "Regge trajectories" of defect local 
operators in terms of the approximate higher spin symmetry expected to hold in the 3d 
Ising model. 

3.1 The set-up 

Let us now describe in some more detail the realization of monodromy defects in the 3d 
critical Ising model on a cubic lattice with periodic boundary conditions. 
The partition function of the model on a cube of side i is 



Z=^exp[-/7(K})] , 
where the sum runs over all the spin configurations. The Hamiltonian reads 



(3.5) 



^(Kl) 



-/3c 2^ Jl^xy) ^x(yy , 



(3.6) 



where the sum runs over nearest-neighbor sites and the Z2 variables Gx are defined on the 
sites. 

The coupling /3c is set to the best known critical value 0.22165455 [15]. In the following 
we shall compare our results with the existing estimates for the spin and energy critical 
dimensions. The most precise estimates for these quantities are Aq- = 0.51813(5) and A^ = 
1.41275(25) from Monte Carlo simulations [16] and A^ = 0.51819(7) and A^ = 1.4130(4) 
from Strong Coupling Expansions [17]. A conservative combination of these results gives 
the two values A^- = 0.5182(2) and A^ = 1.4130(5) which we shall use in the following. 

As anticipated in Section 3, the position of the monodromy defects is encoded in the 
sign of the couplings of the spin-spin interaction J{xy) ■ We set these couplings everywhere 



to +1, except on the bonds that intersect a surface S joining two defect hnes on the dual 
lattice as depicted in Fig. 1, for which we choose J{xy) = ~1- 

On a finite lattice it is not possible to define a single straight defect; the simplest 
choice is to put a couple of defects as far as possible from each other and measure the 
observables in the neighborhood of one of them; if the size of the lattice is large enough the 
distant defect lines do not disturb the measure. In some cases, in the correlation functions 
involving local operators on the bulk and on the defect the presence of other defect lines, 
including the copies generated by the periodic boundary conditions, cannot be neglected. 

Our aim is to realize on the lattice the lowest dimensional operators o living on the 
monodromy defect (this will be the focus of the next section) and to compute their anoma- 
lous dimensions by means of Monte Carlo simulations. We will also study some correlators 
of local operators in the bulk and on the defect, as discussed in Eq. (2.3) 

As a basic update algorithm we chose the standard Metropolis algorithm with multi- 
spin coding technique. Our version of this method is able to update 64 independent lattices 
in parallel on a simple desktop machine. It is important to notice that the defect breaks 
the translational symmetry in the two transverse directions and it would be a waste of 
CPU time to update the whole lattice before every measure. A simple way to speed up the 
simulation is to update sub-lattices of decreasing transverse dimensions centered around 
the defect in a hierarchical way [18]. In order to avoid finite size effects the lattice size 
i, in all the directions, has to be large enough; it turns out that £ = 70 is adequate for 
computations involving even spin operators, and £ = 120 is adequate for odd ones. 

3.2 Representations of the lattice symmetry group 

As discussed in section 3, in a continuum theory with a monodromy line, defect local 
operators have definite scale dimension and S0(2) spin. We will realize some of such 
operators on the lattice in terms of spin operators sitting close to the monodromy line. 
These realizations are classified, rather than by their S0(2) spin, by the representation of 
the corresponding discrete symmetry on the lattice. 

Let us thus consider the symmetry of the lattice plane orthogonal to the defect line, 
as in Fig.s 1,2. In absence of the monodromy cut corresponding to the projection of the 
defect S, the symmetry of the square lattice would be given by the D4 dihedral group, 
generated by the rotation of 7r/2 (say, counter-clockwise) and the reflection about any 
of the symmetry axes (say, the horizontal one in Fig. 2). In presence of the cut, the 
counter-clockwise rotation of tt/2 must be accompanied by a "gauge" transformation to 
bring the defect to its original position; in doing so, the spin variables which switch side 
w.r.t. the frustrated plane change sign; we denote this transformation as a. Its action on 
the elementary spins cjj {i = 1, 2, 3, 4) at the corners of a plaquette "linked" with the defect 
(as depicted in Fig. 2) is as follows: 

a: (cTi, 0-2, 0-3, 0-4) 1-^ (0-2, 0-3, 0-4, -0-1 ) • (3.7) 

Indeed, bringing the cut to its original position after the rotation, it crosses the ui spin 
variable, switching its sign. 




Figure 2: The 2d lattice of a plane transverse to the defect line. The projection of the defect plane 
is the heavier line, crossed by frustrated links, which plays the role of a Z2 monodromy cut. 



The reflection b with respect to the axis through the origin containing the projection 
of the frustrated plane (the x^ axis) is not affected by the presence of the cut and acts on 
the spins o"j as follows: 



b : (cri,Cr2,Cr3,Cr4) l-> (o-4,0-3,C72,C7i' 



(3.8) 



These transformations satisfy a® = 1, 6^ = 1 and {ab)^ = 1, and generate thus the 
dihedral group Dg. Thus, in our discrete lattice set-up, the symmetry group in the plane 
is effectively augmented in presence of the monodromy, and mixes the space-time and the 
flavor Z2 symmetry; Coleman- Mandula theorem does not apply in this case. 

Besides the -Dg invariance, there is another symmetry of the theory which turns out 
to be useful in the classification of the local operators on the defect, namely the reflection 
with respect to a plane orthogonal to the defect line; if x^ is the coordinate along the 
defect, the reflection with respect to such a plane through the origin is of course 



x^ t-^ 



-X 



(3.9) 



and we call the parity of an operator with respect this symmetry 5-parity. 

The Dq group has order 16 and possesses seven irrepses, four of dimension 1 and 
three of dimension 2. The 4-dimensional representation acting on the spins Uj of Fig. 2 
decomposes into two bi-dimensional representations, which we denote as -ffi/2 aiid H^^2- 
A basis for i/1/2 is given by {ijj, ip*), with 



V' 



0"! + a;o"2 + u; as + u (T4 



where uj = exp(i7r/4). In this representation we have 



.r. 




.r. 




(3.10) 



(3.11) 



With our conventions, a field ^p transforming under a as o(/? = e~ '^(/9 lias spin J, so 
i> (ip*) carries spin J = 1/2 (J = -1/2). 

The basis for the representation H^/2 is instead given by (V'3/2) ^'3/2)5 with 

V'3/2 = CTl + UJ CJ2 + OJ CJ3 + aJ(T4 . (3.12) 

In this representation we have 

« fe'=^ = ("" "^ M . o (ly^] = I :, : 1 1 :;.r 1 . (3.13) 








Both the representations -ffi/2 a-^id .^3/2 are odd under the the flavor group Z2. All 
other irreducible representations of Dg are Z2 even and can be obtained by decomposing 
the direct product of iifi/2's. Three of them can be realized in terms of the bilinears UjUj 
corresponding to the four links (ij) of Fig. 2 (see also Fig. 6). Indeed the 4-dimensional 
representation £4 acting on the four links can be decomposed as the sum of a bi-dimensional 
representation, which we call V, and two unidimensional representations: £4 = V(BS(BT^- 
Here S is the trivial representation, which acts on 

s = cri(T2 + 0-20-3 + 1730-4 - 0-4 0-1 (3-14) 

hy as = s,bs = s. The basis elements (-Di, D2) of the two-dimensional vectorial represen- 
tation V, corresponding to spin J = 1, can be chosen to be 

Di = a\a2 + 0-20-3 — 0-30-4 + 0-40-1 , D2 = -0-10-2 + 0-20-3 + 0-30-4 + 0-40-1 . (3.15) 

In this representation the generators act as follows: 

\ / n 1 \ / 7~i \ 

(3.16) 

Notice that the generator a acts on the the combination D = Di + 1D2 as aD = e~"^'^ D, 
so that D has spin 1, while, of course, D = Di — \D2 has spin —1. Finally, T4. is a 
representation of spin J = 2 acting on 

t+ = 0-10-2 - 0-20-3 + 0-30-4 + 0-4 0-1, (3-17) 

by at-i- = — 1+, 6t+ = t+ 

The D^ representations constructed up to now are defined on a single plane orthogonal 
to the defect line, thus are all even under 5-parity (3.9). There are two more unidimensional 
representations of Dg, which we shall denote as P and T_. P is generated by antisymmetric 
products of two representations of -ffi/2 type, and cannot thus be realized in terms of the 
links aiGj of the single plaquette in Fig 2. Their minimal lattice realization involves the 
cube depicted in Fig. 3, obtained by adjoining to the square of Fig. 2 its translation of 
one lattice spacing along the defect line L. We denote the spins of this new square with 
(t'^. Note that the 5-reflection defined in (3.9) now is 

S -.ai^a'i. (3.18) 
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Figure 3: The spin variables around the nionodromy Une in terms of which the defect operators 
of lowest dimensions of table 2 can be realized, as described in the text. 
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Table 1: Schematic description of several lattice operators, built as bilinears in -ip and "03/2 in 
analogy to the free field approximation to the primary operators of the continuum theory, df 
denotes the finite difference df{x) = f{x+l) — f{x) and g d f = gdf—fdg. Using the transformation 
properties of tp and V'3/2 one can verify at once the transformation properties of these bilinears, 
in accordance with the decomposition of representations described in eq. (3.27) and Fig. 5. As 
described in the text, some of these operators can be built explicitly from the spins at the vertices of 
a single plaquette, some require us to use spins from the vertices of a cube. These lattice operators 
provide natural candidates for the corresponding operators in the continuum theory, up to some 
ambiguity due to the fact that the spin J on the lattice is defined modulo 4. 

The anti-symmetric combinations P4 of the diagonals of the four faces which do not 
intersect the defect line define a four-dimensional reducible representation of Dg which can 
be decomposed as P4 = P © F © r_, where now the representations P, V and T_ act on 
5-odd operators, that we call respectively p°,D°, and t?_. We have 



P" = klO-g] + [cr2Cr3] + ['73Cr4] - [a^a^] , 

D'l = -[aia'2] + [a2<T'r^] + [a3<] + [aia[] 

D°2 = [(Tia'2] + [as^] - [^3<] + [f^4^'i] , 

tl = [aia'2] - [ct24] + [asa'^] + [aia[] , 



(3.19) 
(3.20) 
(3.21) 
(3.22) 



The transformation of D° under a and h is the same of that of 



with [aia'j] = a^a'^ — (j[aj. 
D, defined in (3.16), while for the other two operators we have ap° = p°, bp° = —p° and 
at°_ = —t°_, bt°_ = —t°_. Similarly the anti-symmetric combinations P2 of the principal 
diagonals of this cube can be decomposed in the sum D2 = ?+ © P. The pseudoscalar 
representation T^ now acts on 



t\ = [criCTg] - [a2(j'^] 



(3.23) 
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Figure 4: An alternative set of spin variables in terms of which we construct some of the defect 
operators, as described in the text. These spins are lying around the monodromy line in plane 
orthogonal to it, just as in Fig. 2. In particular, the 5-even pseudo-scalar representation P can 
be realized in terms of the bilinears CTiCr.i+i and i7i(Ji+i, indicated in the drawing by the two sets of 
diagonal segments. 

as in (3.17). The pseudoscalar representation P, instead, acts now on the 5-odd operator 



V 



[cTlO-g] -F [(720-^] 



(3.24) 



which is however less efficient in numerical simulations than p° introduced in eq. (3.19). 

It is not difficult to convince oneself that for any irreducible representation of Dg one 
can construct primary operators of both 5-parities. Those described so far can be also 
expressed as bilinears in ^ and/or -03/2 as shown in tab. 1. The other primary operators of 
opposite 5-parity can be expressed as multi-linear products of Oi or as operators involving 
more nodes of the lattice; they are thus difficult to deal with in numerical simulations and 
are expected to have larger anomalous dimensions. 

As an example we describe a simple realization of an 5-even pseudoscalar p which 
can be obtained by considering, instead of the eight spin variables on the cube of Fig. 3, 
eight spin variables lying around the line defect as in Fig. 4. The generators act on these 
variables as follows: 



a: (cri,(T2,(T3,(T4) I-)- (cr2,o-3,(T4,-cri) , (5-i,a2,a-3,o-4) I-)- (0-2,0-3,5-4,-ci-i) ; 

6 : (71 O 0-4 , (T4 O (Ti , (T2 ^ <73 , 0-3 O CJ2 . 

It is easy to check that the Z2-even operator 

p = 0\a2 + Cr2(T3 -|- (73(74 — (7i(74 — {a\a2 + Cr20"3 -|- 0"3(74 — a\a^ 



(3.25) 



(3.26) 



transforms in the pseudoscalar representation P, i.e., we have ap = p and bp = —p. 

An useful tool to summarize the Dq irrepses we discussed above is the graph associated 
with the decomposition of the tensor product of any irreducible representation Ri with the 
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Figure 5: The irrepses of the D8 dihedral group are encoded in the extended Dynkin diagram of 
the Vq algebra, as described in the text. Open circles denote one-dimensional representations while 
grey circles are associated to two-dimensional representations. 



two-dimensional representation H1/2, see Fig. 5. The incidence matrix of this graph^, 
which turns out to correspond to the extended Dynkin diagram of the Pg Lie Algebra, is 
given by the Clebsh-Gordan coefficients Cij in the decomposition 



Hi/2 <E) Ri 



2_^CijRj . 

j 



(3.27) 
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T+,T^ 
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- 
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Table 2: The lowest anomalous dimensions of the local operators at the defect line. The round 
brackets indicate the statistical error, while the square brackets for A > 3 denote an estimate of 
the systematic error (see more details in the text). For the operator f^, realized on the diagonal of 
elementary cubes, we did not use a direct Monte Carlo evaluation which is too noisy, but a lower 
bound obtained by considering only the contributions of the spin-spin two-point functions involved. 



4. Results 

We have seen in section 3.2 that most of Dg representations can be realized using the links 
and the nodes of a plaquette topologically linked with the defect as shown in fig. 2. 

In a first series of numerical experiments we evaluated the correlation function along 
the defect line, between a pair of links of this type located at a mutual distance x. We 
performed an independent simulation for every value of x and for every orientation of the 



^This is very similar to the McKay correspondence between discrete SU(2) subgroups and extended 
Dynkin diagrams of ADE type. 
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pair of links. We then arranged these correlation functions in irreducible representations of 
Ds, according to the prescriptions of Eq.s (3.14, 3.15, 3.17) for the operators s, D and t+, 
in order to extract the exponent of their power- like decay. We expect, for the correlators 
on the defect of such operators, the behavior 

{s{0)s{x)) = ass x-^'^^+constant , {D{0)D{x)) = aoD x~^^^ , (t+(0)t+(a;)) = au x~'^^' . 

(4.1) 
We checked in all cases this power behavior, even if the evaluation of the exponents A 
cannot be very accurate, due to the fact it turns out that A > 2. This implies that the 
correlation function falls off very rapidly and after few lattice spacings the signal is drowned 
in noise (this tendency is already evident for the correlator of the operator D, that has 
Afl = 2, depicted in Fig. 7, even if in this case the fit is still very good). When the 
exponents are A > 3 the accuracy of their determination becomes problematic because the 
value depends on the way we fit the data. This means that these evaluations are affected 
by a systematic error. In order to get an idea of the size of this error, we fitted the data 
with two different procedures, namely to a single power law, like in (4.1), or to the binomial 
a 1 1?' + hjx^^ ^^' taking into account also the contribution of the first secondary operator 
which can contribute; we used the difference between these two determinations of A as a 
rough estimate of the systematic error, which is reported in square brackets in Tab. 2. 

At small values of the spin J, the identification between lattice operators in a given 
Dg representations and continuum operators in appropriate 0(2) representations is clear. 
The operators with high spin in the continuum theory are expected to have large confor- 
mal dimension, and thus give subleading contributions to the lattice operators. The one 
exception is the operator p. Although it has no spin from the point of view of Dsj it is 
built from alternating spins in a way which would closely resemble a J = 4 operator of the 
continuum theory. In the free theory, a J = operator with the quantum numbers of p 
would be a descendant of p°. The J = 4 primary would have dimension close to 5. As the 
lattice correlation function of p° and p is very small, the J = 4 contribution to p must be 
very large, and the numerical estimates correspondingly poor. 

In the case of the scalar s we can conveniently extract A^ from the one-point function, 
which is expected to have the functional form 

(s) = —^ — h constant , (4.2) 

where i is the length of the defect line. Combining this one-point function with the two- 
point function defined in (4.1) we can also extract the universal amplitude ratio asj \/\a^\ 
which turns out to be 0.33(1). 

It is worth noting that Kramers- Wannier duality allows to map the link variables of 
the Ising model into the plaquettes of the dual Z2 gauge theory. Precisely we have, for any 
coupling /?, 

Pij = cosh 2/3 - GiOj sinh 2/3 , (4.3) 

where P^ is the plaquette of the dual lattice orthogonal to the link (ij). Thus correlation 
functions between link variables of Fig. 2 can be written as correlators between staples 
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Figure 6: Using Kramers- Wannier duality one can transform the four link variables of Fig. 2 into 
staples which deform the line defect and can be used to build the displacement operator. 
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Figure 7: The correlation function of the displacement operator. The solid curve is the one- 
parameter fit to gdd/x'^. 



deforming the defect line, as indicated in Fig. 6. In particular the vector representation 
V defined in (3.15) allows to build a discretized version of the displacement operator D 
discussed in Section 2. This operator has protected quantum numbers and in particular 
has dimension Ad = 2. This is nicely confirmed by a one-parameter fit of our numerical 
data, as shown in Fig. 7. 

In another set of numerical experiments we evaluated the correlation function between 
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the spin variables associated to the vertices of the plaquette of fig. 2. In terms of these we 
can construct the correlators of the Z2 odd local operators of semi-integer spin J = 1/2 and 
J = 3/2 defined in (3.10) and (3.12). It turns out that A^ is slightly less than 1, so that 
we can follow the signal of the correlator for many lattice spacings. The statistical errors 
are so small that the quality of the data cannot be appreciated in a plot; the results are 
therefore explicitly reported in Table 3. Using these data we can extract a rather precise 
estimate of the anomalous dimensions of this operator, which is given in Table 2 together 
with those of the other lowest-lying operators. 

One emerging feature of the spectrum of anomalous dimensions is that most of the 
estimated values are rather close to the values expected in a free-field theory. This can be 
seen by comparing the dimensions of the operators of table 1 obtained in our simulations 
with the dimensions that the corresponding bilinears would have in a free field theory, 
where one would have [V'] = 1, [V'3/2] = 2 and, of course, [9] = 1. It turns out that these 
free-theory values represent in almost all cases the nearest integers to our numerical results. 



X Gii2{x) 


X 


Gi,2{x) 


2 0.86752(3) 


11 


0.04529(5) 


3 0.46136(4) 


12 


0.03873(5) 


4 0.28125(4) 


13 


0.03341(5) 


5 0.18892(4) 


14 


0.02925(5) 


6 0.13592(4) 


15 


0.02592(4) 


7 0.10287(4) 


16 


0.02301(5) 


8 0.08066(4) 


17 


0.02071(5) 


9 0.06518(4) 


18 


0.01873(5) 


10 0.05383(4) 


19 


0.01692(5) 



Table 3: The values of the spin-1/2 correlation function Gi/2{x) = Kc{ip{x)-ip* {0)) on the defect 
line. 



We also considered mixed correlation functions of local operators in the bulk and on 
the defect as discussed in Eq. (2.3). In particular, we placed at the origin a defect operator 
Oj of spin J and took as bulk operator the spin a associated to a node at a distance |x*| = r 
from the defect line and at a distance {x'^l = Vd'^ + r^ from the origin, see Fig. 8. In this 
case if the critical Ising model is a conformal-invariant theory we expect that 



{a{d,r,<P)oj{0)) = Cje 



-i</>J_ 



(4.4) 



(^2 + d2)A„j ' 

where A„ = 0.5182(2) is the anomalous dimension of a, while (f> is the azimuth angle 
around the defect line, with (p = corresponding to the surface of frustrated links. The 
setup drawn in Fig. 8 corresponds to (f> = n. In our numerical calculations we observed the 
best signal at J = 1/2, with the operator 01/2 corresponding to ip. Notice that Eq. (3.10) 
fixes completely the phase of this mixed correlation function. In fact we have 



arg((c7((i, r, (p)ij{0))) = arg(l + w^) + 



vr 



7 



-vr ■ 



2 



(4.5) 
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Figure 8: Set-up for the mixed correlation function of the scalar spin on the bulk and the Z2- odd 
local operator on the defect. This local operator is built with the spins associated to the corners of 
the plaquette wrapped around the defect line, following the prescriptions of Eq.s (3.10) and (3.12). 



Similarly, if we considered the operator O3/2 = V'3/2 we would have 



arg((a(d,r,</>)^3/2(0))) = arg(l + a;) + 



Svr 
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-vr 



(4.6) 



We plot in Fig. 9 the numerical data of the imaginary part of ((T(d, r, 7r)^(0)) taken at 
fixed longitudinal distance d = 10, as well as its one-parameter fit^ to (4.4). 

Another observable we studied in our simulations is the one-point function of the 
energy operator in presence of a line defect. This can be considered as a mixed bulk- 
defect correlation function in the case in which the defect operator is the identity. In this 
particular case Ai = and Eq. (2.3) gives 



(6(X)1(0)) 



Cf 



„«|Ae 



(4.7) 



where |x*| is the distance from the defect, while C| is a numerical coefficient which "mea- 
sures" the bulk-to-defect pairing of the energy operator, which is normalized in the bulk 



in the standard way, i.e. (e(x)e(O)) 



-2Aj 



, with A, = 1.4130(5). 



In our simulations we used as a probe the link variable ataj, which can be decomposed 
as the sum of the identity and the tower of local energy operators e, e', . • • • If this probe 
is sufficiently far from the defect line only the operators of lowest dimension contributes; 
thus we expect to have, in presence of a defect line, 



((TiO-j(x))defect = (cricrj(a;)l(0)) = Constant + 



X 



jlA, 



-|- higher order terms , 



(4.8) 



^In comparing our numerical results on the torus to eq. (4.4), written on the covering space, we have 
to take into account the fact that in the covering space there are replicas of the defect lines. The replica of 
the J = 1/2 defect operator closest to the bulk operator a is the one that would appear at distance £ in the 
vertical direction from the one drawn in Fig. 8. With respect to it, the bulk spin is at transverse distance 
r' = £ — r and at an angle (f)' = 0; we find that the corresponding contribution of the form (4.4) has to be 
included in the fit, while further replicas give neglegible contributions. 
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Figure 9: The correlation function between the scalar spin on the bulk and a J — ^ defect local 
operator in a cubic lattice of size i^ with £ = 120 with periodicboundary conditions. The solid line 
is a one-parameter fit to (4.4) with d = 10. In this fit we also considered the contribution of the 
nearest copy of the defect line, see the discussion in footnote 2. 



while in the bulk, in absence of the defect line, the correlator of two link variables aiaj and 
akai is 



{aiaj{x)akcri{0)) = constant + 



+ higher order terms . 



The coefficient C| of Eq. (4.7) is thus given by the universal amplitude ratio 

CI "^ 



(4.9) 



(4.10) 



CJ 



ir"^i 






-0.167(4) 0.968(2) 0.61(9) 



Table 4: Universal amplitude ratios defining the bulk-to-defect pairing associated with the energy 
operator e and the scalar a, defined in Eq. (4.7) and in Eq. (4.11). 



The bulk-to-defect pairing of the scalar field a is encoded in the constant Cj of Eq. 
(4.4) which allows to define further universal ratios 






Cj 



Vc, 



(4.11) 



aa^OjOj 
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where c^a and a^joj are the coefhcients of the 2-point functions of the scalar a on the 
bulk and of the operator Oj on the defect line. In Table 4 we report our estimates for 
these universal ratios for J = 1/2 and J = 3/2 as well as that associated with the energy 
operator. 

5. Conclusions 

The first lesson we draw from our numerical results is that the hypothesis of conformal 
invariance of the monodromy defect in the critical 3d Ising model seems well supported. 
From a theoretical point of view, it would be useful to fully spell out the conditions for 
a generic scale invariant defect in a CFT to be conformal, generalizing the known results 
for boundary conditions [8]. It would also be interesting to test this assumption in other 
concrete examples. 

The second lesson is that numerical methods are suitable to derive information about 
the lowest lying operators on the defect, including both their anomalous dimensions and 
various OPE coefficients. Our results suggest it may be interesting to revisit the problem 
of computing numerically OPE coefficients and correlation functions in the bulk theory, 
and compare them with the results of bootstrap. 

The most immediate direction for future inquiry is to test our results against analytic 
and semi-analytic methods such as the conformal bootstrap and the e expansion. Notice 
that the monodromy defect can be defined uniformly in a scalar theory with (j)"^ interaction 
in the whole interval of dimensions 2 < D < A where the bulk theory is expected to have 
an infrared conformal fixed point. It should thus be possible to compute both anomalous 
dimensions and correlation functions in 4 — e and 2 + e dimensions. 

We find interesting that the spin 1/2 operator on the monodromy defect has dimension 
only slightly smaller than 1, and is thus very weakly relevant when integrated on the defect. 
This raises the possibility of a fully perturbative RG flow to some nearby conformal fixed 
point without rotational symmetry. It would be interesting to investigate this fiow. 
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